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Abstract 

In this paper we reconsider recently derived bounds on MeV tau neutrinos, 
taking into account previously unaccounted for effects. We find that, assuming 
that the neutrino life-time is longer than 0(100 sec), the constraint N e ff < 3.6 
rules out v T masses in the range 0.5 (MeV) < m Vr < 35 (MeV) for Majorana 
neutrinos and 0.74 (MeV) < m Ur < 35 (MeV) for Dirac neutrinos. Given that 
the present laboratory bound is 35 MeV, our results lower the present bound to 
0.5 and 0.74 for Majorana and Dirac neutrinos respectively. 
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Despite a considerable experimental effort it is still unknown whether or not neutrinos 
have a nonzero mass. The following upper bounds are known IJ: 

m Ve < 10 eV (la) 

m t , ll < 270KeV (lb) 

m UT <35MeV (lc) 

More stringent bounds can be determined from cosmological considerations. Assuming that 
there is no cosmological constant, the masses of light (m v < 0(100 Me V)) stable neutrinos 
are bounded by the Gerstein-Zeldovich limit || which can be written as 

/ 1Q 10 \2 

m u < 380eV[—. - h 1Q0 ) , (1) 

\tu{yrs) J 

where tjj is the Universe age and hioo = H/ 100 km/sec/Mpc is the dimensionless Hubble 
parameter. By most estimates tu > 12 Gyr and 0.5 < h < 1. Thus the cosmological 
upper bound the neutrino mass (all weakly interacting flavors) is 40 eV. If there is a non 
zero cosmological constant the bound is somewhat less restrictive, m v < 200 eV 0j. These 
bounds are obtained by ensuring that the energy density of relic cosmic neutrinos be less 
than the closure density, and do not apply if neutrinos are unstable with life-times smaller 
than the the age of the universe. Recently, it has been pointed out [4] that nucleosynthesis 
considerations can further constrain the mass of an unstable r neutrino. 

Nucleosynthesis calculations along with data on light element abundances constrains the 
number of effective neutrino species contributing to the cosmic energy density N v , to be less 
than 3.6 (Note that in ref. [4] a slightly stronger bound N v < 3.4 was used). This bound 
is a consequence of the fact that the rate at which the universe cools depends on the total 
number of species contributing to the cosmological energy density which in turn determines 
the light element abundances. Previously, it was thought that the contribution from a heavy 
neutrino species with m > O(MeV) could be neglected since its energy density was assumed 
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to be Boltzmann suppressed. However, after a massive neutrino decouples and becomes 
non-relativistic its' energy density grows relative to that of a massless species. Therefore, if 
the number density of a massive species at freeze out is on the order of the number density 
of a massless species, then the heavy species may have a greater effect on the light element 
abundance than a massless species. 

In this paper we reconsider the bounds derived in ref. [4] taking into account a few effects 
which were neglected in that paper. Our results for the number densities are in reasonable 
agreement with theirs for the case of massive Majorana neutrinos but differ in the case of 
massive Dirac neutrinos. There is a rather large difference between our respective results 
for the effective number of neutrinos species contributed by the heavy v T . We find a more 
stringent bound which allows no window between the experimental bound and the bounds 
derived here. 

In the standard calculations of the relic abundance of a particle species disappearing due 
to annihilation, the following two essential assumptions are made: The particles in ques- 
tion remain in kinetic equilibrium, and Boltzmann statistics are applicable. In this case 
the Boltzmann kinetic equation can be reduced to an ordinary differential equation. If the 
species considered is relativistic or non-relativistic at the time of decoupling, the equation 
can be greatly simplified and relatively accurate analytical solutions can be found. However, 
for a species which is semi-relativistic when it decouples, there is no known analytical ap- 
proximation and one has to integrate the Boltzmann equation numerically. As such, we will 
present some of the details of the calculations. 

If we are concerned with the abundance of species tp then the Boltzmann equation is 
given by (we specialize to the case ip + 2 — > 3 + 4) 




(2) 
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where H = R/R is the Hubble parameter. In this notation 

n^ = g J -j^-^U(E,t), (3) 

fy(E,t) is the phase space density of species if), and g is the number of degrees of freedom 
for if>. C[f] is the collision term given by 
g f nin d 3 p^ 

|M^ 3+ 4| 2 (^/2(l±/ 3 )(l±/4)) 
- | M 3+4 ^ + 2 I 2 (/ 3 /4(l ± U) (1 ± 2))]], (4) 

and 

The amplitude M is summed over initial and final spins. Assuming that the initial and 
final states are in kinetic equilibrium we may write f a = n a exp{-E/T)/n^ Q . Eq.(4) may be 
written as 



dt 



+ 3 ^ n V> = ~9i>92 J dU^dU 2 



o-v M oiie r exp(-(E 1 p + E 2 )/T) K v \ Q % Q 



(n,n 2 -n E ^) 

EQ 1 

n,,,n 2 



(6) 



Following the notation of ref || we have defined 



VMoller = . ^ ^ = [\v^ ~ V 2 \ 2 - \ X V 2 \ 2 ] 1/2 , (7) 

where F is the particle flux. In the non-relativistic limit the Moller velocity reduces to the 
relative velocity. 

It is convenient to introduce the dimensionless variable x = m/T and the relative number 
density r = u/uq, where m is the mass of annihilating particles and uq is conserved in a 
comoving volume^]. In terms of these quantities the Boltzmann equation takes the form: 



dr tiq , \ / 2 2 \ 

^ = ^T(-Mo^)(r -r EQ ) 



2 Often the number density n is normalized to the entropy density s. This choice is convenient because s 
is usually conserved in the comoving volume. However, if a species is annihilating while it has considerable 
energy density, then when it drops out of equilibrium the entropy will not be conserved. 
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It can be shown that || 

< WMoller > z 



a{s-4,m 2 )^Kx{^/T) 



(9) 



Um? 8m A TK 2 (m/T) 

In this equation K\ and K2 are the modified Bessel functions and s = (p^ + P2) 2 ■ The 
squared amplitude integrated over the final particle phase space, for the annihilation of two 
Majorana neutrinos with mass m„ T into fermions with mass mj is given by 



j dU 3 dU 4 I M 

spins 



2 

M 



8nG 2 F j \(C 2 + C\){ml\m) - s/2)) + {C\ - C 2 v )m 2 f (3ml - s/2) 



+^G 2 F ™ 2 (C V + C 2 A ) 
12 eZ 



w 



[s/2 - m UT ) + — (s + 2m)) 



2 w T/ 4 

The analogous expression for the case of annihilating Dirac neutrinos is given by 
1 
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(C- v + C-i)(s-mj-miw 2 /s 2 



wm 



L (C 2 -C 2 A ) (s/2 -ml). 



We have defined w = (s 2 — 4sm/), where m/ is the mass of the final state particle, which for 
our purposes is the electron, or the other neutrinos species which we will take to be massless. 
In the non-relativistic limit the cross section for the annihilation of Dirac neutrinos reduces 
to 
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2 c 2 



l + /3 2 (o + 



(12) 



2 V" ' '" x 2 ' 2(l-z 2 )\ 

where z = nif/M, mf is the decay product mass and (3 is the velocity in the center of mass 
frame. v re i = 2(5 in the rest frame of the plasma. We include this expression because our 
result differs slightly from those stated previously in the literature ||. 

To perform the numerical integration it is necessary to specify the function n Q (T) as well 
as the function T(t). In the simplest case the energy density is dominated by relativistic 



particles in thermal equilibrium T = —HT and n oc T 3 . In what follows we take no = 
0.181T 3 which corresponds to the equilibrium number density of massless fermions with two 
helicity states. The Hubble parameter is given by the expression 

H = ^8np/3m 2 Pl = 1.66^ T 2 /M Pl (13) 

In this simple case all the quantities are defined in terms of the plasma temperature T 
and the numerical integration is straightforward. When there is a mixture of relativistic 
and nonrelativistic (or semirelativistic) species with a conserved number of nonrelativistic 
particles, one must substitute for g* the expression 

g* = g? el (l + p NR /p Rel ) (14) 

where p Rel and p NR are the energy densities of relativistic and nonrelativistic species respec- 
tively. If we neglect the energy exchange between the relativistic and nonrelativistic species 
the rate of the cooling would be the same, T = — HT and we have to solve kinetic equation 
(8) subject to (14). For the choice of normalization n = 0.181T 3 eq.(14) takes the form: 

9* = + 0.554(£>r/T (15) 

with (E) being the average energy per particle for the massive fermions. The massive 
neutrinos will have a thermal distribution as long as they remain in kinetic equilibrium, which 
is true down to temperatures T m 2 MeV. Below this temperature the distribution function 
in the phase space changes from exp(— ^ {m? + p 2 ) /T) to exp(— ^ (m 2 /T 2 + p 2 /T 2 )). This 
change in the distribution function will modify the Boltzmann equation. For large masses it 
is not expected that these changes will have any significant effect. We have calculated affect 
of such changes and found them to be less the two percent for masses above 1 MeV. For 
masses m < 1 MeV, the neutrinos decoupled while they were still relativistic, so their number 
density will not be suppressed relative to a massless species. Thus in the case of m < 1 MeV, 
r — 1, and the massive neutrinos distribution in phase space is given by exp(— p/T). Given 
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this distribution function the average energy is (E) « (m 2 + OAUmT + 3.151T 2 ) 1 / 2 . This 
approximate expression is accurate up to 0.5% and was used in our numerical work for the 
case of light neutrinos. We should also point out that for smaller masses, m = 1 — 5 MeV 
the Boltzmann approximation is not accurate and the helium abundance will decrease by 



5-10% pi]]. However, as we shall see, this will not affect our bounds. 



Exchange of energy between massless and massive particles gives rise to a faster cooling. 
Covariant energy conservation demands that 

T \ gl el + 0A2rf) gl el + 0A2rf v ; 

This is valid if p NR < p ReL and the massive and massless particles are in thermal contact. If 
this inequality did not hold the cooling rate would be different, and eq.(8) would not be valid 
since it was based on the assumption that the massive particles were in kinetic equilibrium. 

We have numerically solved kinetic eq.(8) for gl el = 9 and g™ r given by eq. (15). This 
corresponds to the case of two light neutrinos and one heavy. Note that g™ r depends on the 
unknown function r(T) which is determined from the solution of the kinetic equation. Our 
results for the freeze out number density of the massive Majorana and Dirac tau-neutrinos 
are presented in fig. 1. We then calculated the n/p ratio taking into account that r is not 
a constant but decreases due annihilation of the heavy neutrinos. Furthermore, we used the 
exact expression for the average energy density of the heavy neutrinos by integrating the 
distribution function. 

Bounds on the neutrino mass are derived by computing the net increase in helium pro- 
duction due to the presence of the heavy neutrino species. In accordance with the standard 
nucleosynthesis calculations almost all neutrons which survived down to the temperature 
T 7 = 0.065 MeV turn into He 4 . Thus, by comparing the neutron to proton ratio calculated 
for two massless neutrinos and one neutrino with mass m UT , to the ratio calculated for a 
variable number of massless neutrino species we may bound the neutrino mass. If the n/p 
ratio yield for a given mass m UT exceeds the ratio calculated for 3.6 massless species, then 
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that mass is ruled out. For the case of Dirac neutrinos we assume that the right handed 
species will be populated if tud < 0.74 MeF|12[]. From these considerations we find that, 
assuming the neutrino lifetime is longer than 0(100) sec, the following limits apply 



0.5 < m M < 35 (MeV), (17) 
0.74 < m D < 35 (MeV). (18) 

Tau neutrinos with their mass in the region considered must be unstable in accordance 
with the Gerstein-Zeldovich limit. For our bounds to be valid the neutrino must not decay 
prior to primordial nucleosynthesis. Since the characteristic temperature scale is near 0.1 
MeV the life-time should be larger than or of the order of 100 sec. For life-times shorter 
than 1 secf] , there will be no bound for Majorana neutrinos, but for Dirac neutrinos there 
will be small region of excluded masses near 1 MeV . 

One could deduce bounds on both m VT and on t Vt accounting for the decay of v T in the 
kinetic equation governing both the number density of v T and the n/p-ratio. The final n/p- 
ratio depends not only the life-time of v T but also on the type of particles in the final state. 
A considerable effect associated with the decay might emerge from the distortion of the 
electron neutrino spectrum if there is a decay into v e . This is analogous to the distortion of 
the spectrum of the electron neutrinos due to electron-positron annihilation in the standard 
scenario at the level about 1% found recently PlJlOfl. The size of the effect in the case of 
massive v T annihilation or decay should be bigger because the (hypothetical) mass of v T is 
assumed to be larger than m e = 0.5 MeV. We hope to take all these effects into account in 
the subsequent publication. 
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